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Amendments to the Specification: 

After line 7, page 1 1 of the as-filed application, please insert the following 
paragraphs: 

For example, and as discussed in this reference, multi-variable polynomials have been 
successfully employed to realize functional index fits. It is well known that the quality of the 
interpolation does not increase as the order of the approximation increases. C. Runge 
examined a famous example in one dimension, on the interval -5<x<5 where he showed that 
the maximum of the error approached infinity as the polynomial order increased. 

To avoid such problems, the domain may be partitioned into smaller sub-intervals. A 
piecewise polynomial interpolation is then obtained over these smaller intervals. Perhaps the 
most popular is a cubic; both the first and second derivatives are continuously differentiable. 
When the approximating polynomial g(x) of the original function f(x) satisfies the property 
that g=f and g'=f ' at the end points of the subinterval, then 

b b 



This simply infers that the "strain" energy minimized by the spline approximates well 
the energy of the original function. That fact along with the ease of computing the partial 
derivatives in any one direction make a multi variable spline an ideal choice for an 
approximating function. 

For example, consider a function fit to only one variable. The truncated Taylor series 



provides a satisfactory approximation to f(x) provided that x is sufficiently smooth and f is 
sufficiently close to a. If the interval is large, requiring large values of n, spurious oscillatory 
problems may result. This of course precipitated the use of splines in which the interval was 
truncated and fit with polynomials of lower degree each with several continuous derivatives 
to generate smooth functions. The pp-form of a spline provides a description of the functional 
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in terms of intermediate breakpoints (or knots) £ g m and coefficients c u to be determined 
for a spline of order k so that 

For this analysis, the B-form of the spline is used where 

/(*) = Z^# M (jt). (22) 

j=i 

where Bj k is the j* B-spline of order k with knot sequence t h t 2 t n + k . The B form has 

become the standard way to represent a spline during its construction. The knots refer simply 
to the sub-intervals. If f was related to 3 variables x, y, and z, the trivariate tensor spline 
approximation would be 

/Uy,z) = HI c{q,r,p)B rk {x)B qJc {y)B pk {z). (23) 

q r p 

Note that as additional variables are incorporated, they are handled as a product of the 
polynomials used to work with one dimension. The task is to determine the trivariate 
coefficients c. For this purpose consider f defined as sin(x) * (cos(y) * tan(z)) 2 . The variables 
x, y, and z are allowed to vary from 0 to 2nl5 each with 20 intermediate points in a nested do 
loop. The fit of these 8,000 points with only 729 coefficients has a max error of 0.0291. This 
function is smoothly differentiable. 



3 



DOCKET NO. NNI-0052 

Application No. 09/873,622 

Office Action Dated: January 24, 2006 



PATENT 



After line 12, page 11 of the as-filed application, please insert the following 
paragraphs: 

For example, and as described in this reference, a performance objective may be to 
choose b (thickness of a backiron) and t (thickness of a magnet) to maximize the average 
magnetic field in a 4" air gap divided by the weight of a magnetic vehicle. The magnets have, 
for example, a residual B field of 1.25T and a working coercivity of 955200 A/m. The carbon 
steel backiron thickness b may range from 0.5 - 10 inches in thickness, in a total of ten steps. 
The magnet thickness t may be stepped through the range of 1 - 12 inches with a step size of 
one inch. The total number of test cases to cover the range of values dictated by both of these 
quantities changing at the same time is 120. 

With an initial backiron thickness of 0.5 inches, everything above the interfacial line 
separating the magnet and the iron may be scaled vertically from a range of 1 - 20. 
Everything below the interfacial line between the magnets and the iron may be scaled from a 
range of 1 - 12. To maintain an air gap of 4", the mid point of the air gap may also be scaled 
in step with the magnet to maintain the separation distance of 2 M . Thus, there are a total of 
three parametric scaling loops for this analysis, one for the backiron, one for the magnets, and 
one to adjust the horizontal Neumann Boudnard representing the middle of the air gap. 

We seek the solution of the vector potential due to source currents using fictitious 
surface currents placed on the interface separating dissimilar materials. The equation 
dictating the magnetic field involves a single integral equation for the unknown fictitious 
surface current, 



where jj, = — -andB t is the tangential B field due to all other external magnetic sources. 
Mi 

Equation (24) must be altered slightly to account for media saturation. Additional 
volumetric sources equal to VxM must in general be distributed throughout the volume of 
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the iron. Their magnitude is smaller than the surface currents and is indeed zero when the 
media is linear. The general boundary element equation used for saturable media is 

B = B 0 +$V x(^(r')G(r ) r')) ds'+fv x(J m (r' ) G(r, r' )) dV (25) 

r y 

Obviously these volume currents, J m = V xM , are unknown a priori. Solution of the 
problem is as follows. First, assume J m = 0; second, compute the B field everywhere; third, 
check all saturable regions to determine whether VxM ^0; and substitute J m - VxM. 
The process is then repeated until no change in these magnetization currents is witnessed. 
Under relaxation can be utilized to accelerate conversion, with an under relaxation parameter 
of about 0.3. Throughout this process the magnetization current is always a source term and 
is treated as known. It is computed after the fact and forms a contribution to the right hand 
side of the equations. 

Once computed, the two dimensional variation of median air gap field with magnet 
thickness and backiron thickness must be accurately modeled. A fairly good solution was 
found by first considering the case where the backiron thickness was equal to 0. The working 
magnetic field in the air gap for this case can be modeled quite accurately using a quadratic 
polynomial fit as 

3 

B c Vg = Z C i* 3_ ' = C l' 2 + C 2 ( + C 3 (26) 

Now the issue arises as to how the introduction of magnet backiron begins to 
complicate the determination of the working B field. Of course, the expected magnetic field 
in the gap will increase with backiron thickness. However, the advantage gained by adding 
backiron is somewhat inversely proportional to the magnet thickness itself, since with very 
thick magnets, the iron begins to saturate and becomes less effective for carrying flux. It was 
thus reasoned that any adaptation of equation (26) must be realized using a parameter which 
reflects this effect, one such parameter being simply the backiron thickness "t" divided by the 
magnet thickness "b". A secondary polynomial fit was performed in an attempt to account for 
the enhancement in the magnetic field expected by the addition of backiron as a multiplier on 
the formula found when no backiron exists i.e., 

5 



DOCKET NO. NNI-0052 

Application No. 09/873,622 

Office Action Dated: January 24, 2006 



PATENT 



1=1 



' 'bY (b 

\ 



+ c 5 - +c 6 



(27) 



A least square polynomial analysis was used to determine the parameters as 
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The polynomial match is usually good to less than 5%. 
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